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ABSTRACT 

We numerically investigate the excitation and temporal evolution of oscillations in 
a two-dimensional coronal arcade by including the three-dimensional propagation of 
perturbations. The time evolution of impulsively generated perturbations is studied 
by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-/3 ap- 
proximation. As we neglect gas pressure the slow mode is absent and therefore only 
coupled MHD fast and Alfven modes remain. Two types of numerical experiments are 
performed. First, the resonant wave energy transfer between a fast normal mode of the 
system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, 
the fast wave with global character transfers its energy to Alfvenic oscillations localized 
around a particular magnetic surface within the arcade, thus producing the damping 
of the initial fast MHD mode. Second, the time evolution of a localized impulsive ex- 
citation, trying to mimic a nearby coronal disturbance, is considered. In this case, the 
generated fast wavefront leaves its energy on several magnetic surfaces within the ar- 
cade. The system is therefore able to trap energy in the form of Alfvenic oscillations, 
even in the absence of a density enhancement such as that of a coronal loop. These 
local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of 
wave energy trapped by the system via wave energy conversion strongly depends on the 
wavelength of perturbations in the perpendicular direction, but is almost independent 
from the ratio of the magnetic to density scale heights. 

Subject headings: MHD - Sun: corona, Sun: oscillations, Sun: magnetic fields - waves 



1. Introduction 



The presence of waves and oscillations in the solar corona is a well known feature that has been 



observed for long time. For an overview of the early observational background see Tsubaki (1988). 



Nowadays, because of the increasing spatial and temporal resolution of the EUV instruments on- 
board TRACE, SOHO and HINODE spacecraft, accurate observations of oscillations in different 
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coronal structures are accomplished. Many authors have reported observations of transversal coro- 



nal loop oscillations from both ground and space-based instruments (Aschwanden et al. 1999 



Nakariakov et al. 


1999 


Aschwanden et al. 


2002 


Schrijver et al. 


2002 


). When these observations 


are compared with theoretical models ( 


Roberts et al. 


1984; 


Nakariakov et al. 


1999, 


Nakariakov & 



Ofman 2001), the possibility of inferring some plasma parameters, otherwise difficult to measure, 



and of improving the existing theoretical models is open; see Banerjee et al. (2007) for a review. 
Magnetohydrodynamics (MHD) is the underlying theory of coronal seismology and it is believed 
that all these observed oscillations and waves can be interpreted theoretically in terms of MHD 
modes of different coronal plasma structures. 

The theoretical study of these oscillations and waves can be done from several points of view. 
The first approach is to make a normal mode analysis of the linearized MHD equations, which allows 
to obtain the spatial distribution of the eigenmodes of the structure together with the dispersion 
relation oj(k). Once the elementary building blocks of the MHD normal mode theory are described, 
the main properties of the resulting MHD waves can be outlined. Many authors have explored 
the normal modes of coronal structures, beginning with very simple cases such as the straight and 



infinite cylinder 


(Edwin & Roberts 


1983). 


Goossens et al. 


( 


1985 


); 


Poedts et al. 


(1985) 


spectrum of ideal MHD. 


Oliver et al. 


(1993 



of modes in potential and nonpotential arcades. More complex configurations, such as sheared 
magnetic arcades in the zero-/3 plasma limit, have been studied by Arregui et al. ( 2004a|b ) . Other 



authors have studied eigenmodes in curved configurations with density enhancements that represent 



coronal loops (e.g., Van Doorsselaere et al. 


2004, 


Terradas et al. 


2006, 


Verwichte et al. 


2006a|b|c 


Diaz et al.||2006 Van Doorsselaere et al. 


2009 


)• 



An alternative approach is to obtain the time dependent solution of the MHD equations. Using 
this method, Cadez Sz Ballester ( 1995a|b ) studied analytically the propagation of fast waves in 
a two-dimensional coronal arcade for a particular equilibrium, namely one with uniform Alfven 
speed. Oliver et al. (1998) studied the effect of impulsively generated fast waves in the same 



coronal structure. Del Zanna et al. (2005) studied the properties of Alfven waves in an arcade 
configuration, including the transition region between the photosphere and the corona. Other 
studies have analyzed the effect of the loop structure on the properties of fast and slow waves in 
two-dimensional curved configurations (see, e.g., Murawski et al.|[2005 Brady Arber|2005 Brady 



et al. 2006; Selwa et al. 2006 2007), see Terradas (2009) for a review. 



The main aim of this paper is to analyze the effect of including three-dimensional propagation 
on the resulting MHD waves as a first step before considering more realistic situations like the one 



observed by Verwichte et al. (2004), where the effect of three-dimensional propagation is clear. In 
our model there is no density enhancement like that of a loop and the zero-/3 approximation is 
assumed, so only the fast and Alfven modes are present. We focus our attention on the mixed 
properties displayed by the generated MHD waves that arise due to the coupling when longitudinal 
propagation is allowed. The paper is arranged as follows. In §[2] we briefly describe the equilibrium 
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configuration as well as some of the approximations made in this work. In §[3] we present our deriva- 
tion of the linear ideal MHD wave equations with three-dimensional propagation of perturbations. 
In § [4] the numerical code used in our study is described, together with several checks that have 
been performed by solving problems with known analytical or simple numerical solution. Our main 
results are shown in § [5j where the linear wave propagation properties of coupled fast and Alfven 
waves in a two-dimensional coronal arcade, allowing three-dimensional propagation, are described. 
Finally, in § [6] the conclusions are drawn. 



2. Equilibrium configuration 

We model a solar coronal arcade by means of a two-dimensional potential configuration con- 



tained in the xz-plane in a Cartesian system of coordinates (see Oliver et al. 1993). For this 
y-invariant configuration the flux function is 

A{x,z) = BA B cos (j^-) exp (~J^\ (!) 
and the magnetic field components are given by 



z 



B x (x, z) = Bcos ex P i \ 



x \ i 



B z (x,z) = -Bsin f — I exp f -— 1 . (2) 

In these expressions A# is the magnetic scale height, which is related to the lateral extent of the 
arcade, L, by A# = 2L/tt, and B represents the magnetic field strength at the photospheric level 
(z = 0). The overall shape of the arcade is shown in Figure [TJ 

In this paper gravity is neglected and the /3 = approximation is used for simplicity. Therefore, 
the equilibrium density can be chosen arbitrarily. We adopt the following one-dimensional profile 



Po{z) = Poexp 



(3) 



where A is the density scale height and po is the density at the base of the corona. As shown by 
Oliver et al. (1993), the combination of magnetic field components given by Equation Q with the 



density profile given by Equation ([3]) leads to a one-dimensional Alfven speed distribution in the 
arcade that can be cast as 



v A (z) = v ao exp 



-(2-5) 



2A 



B 



(4) 



Here 5 = ^fi- represents the ratio of the magnetic scale height to the density scale height and vao is 
the Alfven speed at the base of the corona. The 5 parameter completely determines the behavior 
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of the Alfven speed profile and hence the wave propagation properties. The case 5 = 2 represents a 
uniform Alfven speed model, while 5 = corresponds to an exponentially decreasing Alfven speed 
in a uniform density configuration. Other values of 5 represent situations in which both the Alfven 
speed and density depend on height in a different manner. 



3. Linear waves 



In order to study small amplitude oscillations in our potential arcade the previous equilibrium 
is perturbed. For linear and adiabatic MHD perturbations in the zero-/3 approximation the relevant 
equations are 



Po^ = —(V x B) x Bi + —(V x Bi) x B, 
dt fi Po 



(5) 



0Bi 

~dT 



V x (vi x B), 



(6) 



where no is the magnetic permeability of free space and the subscript "1" is used to represent 
perturbed quantities. These equations are next particularized to our two-dimensional potential 
arcade equilibrium. As the equilibrium is invariant in the y-direction, we can Fourier analyze all 
perturbed quantities in the y-direction by making them proportional to exp(ik y y). In this way, 
three-dimensional propagation is allowed and each Fourier component can be studied separately. As 
a result of this Fourier analysis the perpendicular perturbed velocity and magnetic field components 
appear accompanied by the purely imaginary number i = \J—1. This is undesirable from a practical 
point of view, since Equations ([5j and ^ will be solved numerically and the code is designed 
to handle real quantities only. Nevertheless, by making the appropriate redefinitions, namely 
vi y = iviy and B\ y = iBi y , it turns out that our wave equations can be cast in the following form 



dt 

~dT 

dviz 
dt 

dB lx 

dt 

dB ly 

dt 
dB u 

dt 



1 



1 

Mopo 

1 



dB lx dB 
dz 



lz 



dx 



B : 



dB 



dB 



>.r . , B z 

dx dz 



la 



+ k y (Bi x B x + B\ Z B Z ) 



dB lx dB u 
dz 



dx 



B x 



d 



-k y vi y B x - — {v lz B x - v Xx B z ) , 

d d 
q~ z (.viyB z ) + — (v ly B x ) , 

d 

-k y vi y B z + — (vi z B x - v lx B z ) . 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



- 5 - 



These equations constitute a set of coupled partial differential equations with non-constant 
coefficients that describe the propagation of fast and Alfven waves. As the plasma (3 = 0, slow 
waves are excluded from the analysis. When k y = 0, equations (p7j) — (|12[) constitute two independent 
sets of equations. The two equations for v\ y and B\ y are associated to Alfven wave propagation. 
On the other hand, the four equations for the remaining variables, v\ x , v\ z , B\ x , B\ z , describe the 
fast wave propagation. The basic normal mode properties of fast and Alfven modes in a potential 
arcade with k y = are described in Oliver et al. (1993), while the case k y / has been considered 



by Arregui et al. (2004a). The time dependent propagation for k y = was analyzed by Terradas 



et al. (2008). When longitudinal propagation of perturbations is allowed (k y ^ 0), the six equations 



and their solutions are coupled so we may anticipate fast and Alfven wave propagation to display 
a mixed nature, in an analogous way to the mixed character of eigenmodes obtained by Arregui 



et al. (2004a) in their analysis of the normal modes of the present equilibrium with k y ^ 0. In the 



following the tildes in v\ y and B\ y are dropped. 



4. Numerical method and test cases 

4.1. Numerical method 

The set of differential equations (|7|)-(|12[) is too complicated to have analytical or simple nu- 
merical solutions except for simplified configurations and under particular assumptions. For this 
reason we solve them by using a numerical code, although comparisons with known wave properties 
have been carried out whenever possible. 

When considering a potential arcade as the equilibrium magnetic field, it is advantageous to 
use field-related components instead of Cartesian components in order to characterize the directions 
of interest related to the polarization of each wave type. The unit vectors in the directions normal, 
perpendicular, and parallel to the equilibrium magnetic field are given by 

VA B 

e ™ = p^Tj' e ± = e y , e|| = — , (13) 

where A is the flux function given in Equation ([I]). These unit vectors are related to the Cartesian 
ones as follows 

(B z e x -B x e z ) „ „ „ (B x e x + B z e z ) 
Gn = MB] ' e± = 6y ' 6|! = j~B~j ' ^ ' 

with | B |= (B x + B^) 1 ! 2 . In the absence of longitudinal propagation (i.e. for k y = 0), these three 
directions are associated with the three types of waves that can be excited, namely v\ n for fast 
waves, vix for Alfven waves, and for slow waves. 
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Since we want to model a coronal disturbance with a localized spatial distribution we have 
considered as the initial condition a two-dimensional Gaussian profile given by 



v\ = v s exp 



(x 



X s 



+ (* 

r,2 



Zsf 



(15) 



where v s is the amplitude of the velocity perturbation, x s and z s are the coordinates of the per- 
turbation's center, and a is the width of the Gaussian profile at half height. In the following we 
use v\ = v\ n to excite fast waves and x>\ = v\ y to excite Alfven waves. When k y = the fast mode 
produces plasma motions purely normal to the magnetic field, while the Alfven mode is charac- 
terized by a purely perpendicular velocity component. When propagation along the y-direction is 
considered, pure fast or Alfven modes do not exist and both produce motions in the normal velocity 



component as well as in the perpendicular velocity component (Arregui et al. 2004a). It must be 



noted that the numerical code solves the time-dependent equations in Cartesian coordinates and 



so the solution has to be transformed following expressions (14) to the field- related coordinates. 



The same applies to the initial perturbation, which must be transformed into the corresponding 
Cartesian components. 



The numerical code (see Bona et al. 2009, for details about the method) uses the so-called 



method of lines for the discretization of the variables and the time and space variables are treated 
separately. For the temporal part, a fourth-order Runge-Kutta method is used. For the space 
discretization a finite-difference method with a fourth-order centered stencil is choosen. For a given 
spatial resolution, the time step is selected so as to satisfy the Courant condition. As for the 
boundary conditions, as we computed the time evolution of two initial perturbations, two kinds 
of boundary conditions are used. First, when the initial perturbation in v\ n is the fundamental 
normal mode of the k y = problem, for the v\ n component line-tying conditions are chosen at all 
boundaries, while for the v\ y component, flow-through conditions are selected except at z = where 



line-tying condition is used. On the other hand, when an initial perturbation like ( 15 ) is considered, 
the large photospheric inertia is accomplished by imposing line-tying boundary conditions at z = 0. 
In all other boundaries flow-through conditions are used so that perturbations are free to leave the 
system. In order to increase numerical stability, fourth-order artificial dissipation terms are included 
in the numerical scheme. In all the simulations the effects of this artificial dissipation have been 
checked to ensure that they do not affect the obtained solution, but just contribute to eliminate 
undesired high-frequency numerical modes. 



4.2. Test cases 

Some preliminary tests have been performed in order to figure out the appropriate values of 
numerical parameters, such as the grid resolution or the numerical dissipation, on the obtained 
results for fast and Alfven waves. The first test we have conducted has been to run the code with 
no perturbation at all and to check that the structure remains stable. The results of this numerical 
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run were completely satisfactory. Then the propagation of linear fast and Alfven MHD waves in a 
potential coronal arcade has been considered. 



4-2.1. Fast wave 



The temporal evolution of impulsively generated perturbations with rather similar conditions 



has been accomplished by several authors: Cadez Sz Ballester (1995a) obtained analytical expres- 



sions for the temporal evolution of perturbations when a coronal arcade is taken as the equilibrium 
state; 



Cadez & Ballester 



(1995b) 



Oliver et al. 



(1998) numerically computed such solutions when 



different initial perturbations are used; and more recently Terradas et al. (2008) showed the main 



properties of the time evolution of fast and Alfven waves in low-/3 environments. These works 
facilitate the comparison of our numerical results with known results as well as with analytical 



ones. 



As shown by Terradas et al. (2008), when different resolutions are used time dependent results 



reveal that the grid resolution in our two-dimensional domain is not a critical factor for the proper 
computation of fast waves and that a good representation of the temporal evolution of perturbations 
can be achieved even with a rather modest resolution of 50 x 50 grid points in the (x, z)-plane. 
As mentioned above, numerical dissipation is introduced in our code in order to ensure numerical 
stability. This dissipation is proportional to an adjustable parameter, or dissipation factor, a n . We 
have conducted numerical simulations for different values of the dissipation factor and it turns out 
that the temporal evolution of fast wave perturbations is not modified. 



4-2.2. Alfven wave 

The properties of Alfven continuum normal modes in a potential coronal arcade described 



by Oliver et al. (1993) allow us to anticipate and identify possible sources of difficulties in the 
numerical computation of Alfven wave solutions. First of all, since they are oscillatory solutions 
strongly confined around given magnetic surfaces (both when propagating or in their standing 
mode version), spatial scales quickly decrease with time and so we can expect a rather important 
dependence of the numerical solutions on the number of grid points used to cover the area in and 
around the excited magnetic surfaces. The situation becomes even worse if we take into account 
that computations in a Cartesian grid do not allow us to locate all the grid points along magnetic 
surfaces. This fact affects the numerical results and adds a numerical damping. Furthermore, when 
time-dependent simulations are considered the sampling rate is no more an independent parameter. 
When the spatial resolution of the grid is defined, the Courant condition gives a maximum value 
for the temporal resolution which in turn sets the maximum frequency that can be resolved. 

We have first generated Alfven waves in our potential arcade model by considering an impulsive 



initial excitation of the v\ y component given by Equation (15) with z s = 1 and x s = 0. This 
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implies that the initial disturbance is even about x = and so odd Alfven modes are not excited. 



As described by Terradas et al. (2008), the spatial resolution of the numerical mesh affects the 



obtained amplitude and frequency values. Better resolution provides a closer value to the analytical 
frequency and less numerical damping. We have also checked the influence of numerical dissipation 
and the results show that only the amplitude, and therefore the damping time, decreases when the 
a n parameter is decreased. The spectral analysis of these oscillations at different heights in the 
structure is shown in Figure [2^,. The resulting power spectrum is compared to the Alfven continuum 



frequencies obtained by Oliver et al. (1993). The frequency associated with the generated Alfven 



waves coincides with the theoretical normal mode frequencies of the system, which gives us further 
confidence on the goodness of our code. Alfven waves stay confined to the vertical range of magnetic 
surfaces that were excited by the initial disturbance, since they cannot propagate energy across 
magnetic surfaces. The initial perturbation is decomposed by the system in a linear combination 
of normal modes, but keeping the even parity of the initial disturbance with respect to x = 0, so 
energy is only found in the fundamental mode, the second harmonic, etc. 

In order to better isolate and show the possible numerical artifacts that the code introduces 
into the numerical solution we have considered a simpler case, the excitation of a particular Alfven 



mode around a magnetic surface. According to Oliver et al. (1993), Alfven normal mode solutions 
can be obtained analytically when 5 = 0. For this reason we now select 5 = 0. 



0. For this reason we now select 5 
The initial excitation could now be given by 



(16) 



vi y (x, z) = v\ y {x)5 A(x, z) - A(x = 0, z m ) 

where v\ y is the regular part of the solution, A(x,z) is the flux function defined by Equation ([TJ), 
and z m gives the maximum height of the magnetic field line in which the normal mode is excited. 
It is important to note that the regular solution has a well-defined parity with respect to the x- 
direction depending on whether n is chosen even or odd. However, since a delta function is difficult 
to handle from a numerical point of view, our normal mode-like excitation is performed by an initial 
perturbation of the form 



v ly \ 



x,z 



vi y (x) exp 



A(x, z) — A(x = 0, z n 



(17) 



For the regular part, vi y (x), the fundamental mode with one maximum along the field lines has been 
chosen. It should be noted that the width, a, of the initial perturbation now causes the excitation 
of several Alfven modes in a set of neighboring magnetic surfaces. It is important to consider an 
initial velocity profile which is sufficiently localized in the direction transverse to magnetic surfaces 
so that only a few of them are excited. As we concentrate on the dynamics of a restricted number 
of field lines around a magnetic surface the consideration of other models, with different values of 5, 
would change quantitatively the generated frequencies, but not the overall qualitative conclusions 
shown here. 



- 9 - 



Figure J2ja shows the temporal evolution of the excited v\ y component at a particular location 
as a function of time for three different values for the width of the initial disturbance. It is clear 
that three different solutions are obtained. The two corresponding to the largest widths are rather 
similar, but the one for the smaller width shows a strong damping. It must be said that the 
exact solution of this ideal system should display no time damping, hence we assert that this is a 
numerical effect, that cannot be attributed to a real physical damping mechanism. This undesired 
effect is less important for larger widths of the initial perturbation since, for a given number of grid 
points, the initial condition is better resolved spatially. 

We next fix the width of the initial disturbance, a, and vary the spatial resolution in our 
domain. Figure [3] shows several numerical simulations when an initial normal mode-like excitation 
(Equation [IT]) is made at different heights. It is clear that larger spatial resolution provides the 
more accurately the undamped oscillatory solution. 

Also from this analysis we conclude that the spatial resolution is not a factor that should 
be taken into account in an isolated manner when considering the numerical description of Alfven 
waves on given magnetic surfaces. Indeed, and because of the Cartesian distribution of grid points in 
a system of curved magnetic field lines, low-lying magnetic lines are poorly resolved when compared 
to high-lying magnetic lines for a given grid resolution. This has implications that are worth to be 
taken into account as can be seen in Figure [3} If we compare signals in Figure [3| we can see that, all 
parameters being the same, closer results to the analytical solution are obtained for higher magnetic 
field lines. We can therefore assert that for the numerical simulation of Alfven wave properties the 
resolution of the grid is an important parameter and that it becomes more critical for low-lying 
magnetic field lines than for higher ones. It should be noted that the conclusions of these tests can 
also be applied to the case in which an impulsive excitation is set as the initial perturbation. 



Numerical results 



In this section we present the main results from our numerical investigation. For simplicity, 
first, the temporal evolution of a normal mode- like fast disturbance is analyzed in order to show 
how and where resonant absorption, due to three-dimensional propagation of perturbations in a 
non-uniform medium, takes place. It turns out that previous results obtained for the normal modes 



of coupled fast and Alfven waves in a potential arcade by Arregui et al. (2004a) can guide us to 



understand the time evolution of the system and the energy transfer between resonantly coupled 
modes. Then, a more complex situation is considered by analyzing the time evolution of the initial 



perturbation given by Equation (15). It should be noted that our first normal mode time evolution 



analysis has been proof very useful to further better understand the resulting coupling process 
between both velocity components when a localized impulsive disturbance is used. 
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5.1. Resonant damping of fast MHD normal modes in a potential arcade 

In order to gain some insight into the propagation properties of coupled fast and Alfven waves 
in our configuration, we first study the time evolution caused by an initial disturbance having the 



spatial structure of a fast normal mode for k y = (propagation in the xz-plane). As shown by Oliver 



et al. (1993), pure fast modes in a potential arcade are characterized by a global spatial structure 
determined by the wavenumbers k x and k z , which give rise to smooth distributions with a given 
number of maxima in the x- and z-directions. This results in a discrete spectrum of frequencies. 



The frequencies and spatial structure of the fast modes with k y ^ were computed by Arregui 



et al. (2004a), who showed that perpendicular propagation produces the coupling of the fast normal 
modes to Alfven continuum solutions, resulting in modes with mixed properties. We have chosen 
as initial perturbation the velocity perturbation v i„, of the fundamental fast mode for k y = 0, with 
one maximum in each direction in the xz-plane. When k y = this produces a standing harmonic 
oscillation of the system, as in an elastic membrane. When k y ^ 0, this initial perturbation is not 
a normal mode of the system, but we expect that, the obtained temporal evolution will not differ 
very much from the actual normal mode of the coupled solution. 

Figure [4] displays the results of such simulation. The first frame for v\ n shows the initial spatial 
distribution of the perturbation. Initially, v\ y , the velocity component associated to Alfven waves, 
is zero. As time evolves, a non zero v\ y component appears because of the coupling introduced by 
the three-dimensional propagation. The panels for v\ y in Figure [4] show that unlike v\ n the excited 
transversal perturbations are not globally distributed in the potential arcade, but only at preferred 
locations, around a few magnetic surfaces. When the v\ n and v\ y signals are measured at one of 
those locations, x = 0, z/L = 0.35, it is seen that the amplitude related to the fast-like perturbation 
decreases in time, while the amplitude of the Alfven-like component of the perturbation increases 
in time, see Figure [5j This is an indication of the wave energy transfer due to the resonant coupling 
of the excited fast normal mode to the Alfvenic solution around the excited magnetic surface. For 
long times a decrease in the amplitude of the velocity component, v\ y , can be appreciated and is 



attributed to numerical damping, for the reasons explained in § 4.2.2 



Further confirmation of the resonant wave energy transfer occurring between the modes can 
be obtained by computing the time evolution of the total energy density in our system. This total 
wave energy can be computed as 



SE(r,t) 



P0 (v 2 lx + v\ y + v 2 lz ) + —{B\ x + B\ y + B\ 



Mo 



(18) 



The right-hand side panels in Figure [4] show the spatial distribution of this quantity as a function of 
time. The different frames clearly indicate that, initially, the energy is distributed globally around 
the center of the system, such as corresponds to the initial perturbation we have used. At later times, 
this energy is transferred to magnetic surfaces around the particular magnetic field line in the arcade 
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where the signals in Figure [5] have been measured. The location of this energy deposition is not 
an arbitrary one. As previous theoretical works on the resonant energy transfer have shown, (e.g., 
Wright|[T992| |Halberstadt Goedblo"edl|T993l |Ruderman eTaI"1[T997l |Arregui et aT1|2004a[ |Russell 



&: Wright 2009), global fast modes resonantly couple to localized Alfven continuum modes at the 



magnetic surfaces where the frequency of the fast mode matches that of the corresponding Alfven 
mode. In our case the spectral analysis of the wave energy densities associated to the normal and 
perpendicular components, plotted in Figure [6j allow us to confirm the resonant energy transfer at 
the location where the fundamental fundamental fast mode frequency crosses the Alfven continuum, 
that exactly corresponds to the magnetic surface where Alfvenic oscillations are excited and energy 
transfer occurs, see Figure [BJd. Although the fast mode frequency crosses other Alfven continua, 
coupling can only occur if the parity of the fast and Alfven eigenfunctions along the field lines is 
the same, see further details in Arregui et al. (2004a). This prevents the coupling with Alfven 
continuum modes with two extrema along field lines. Even if the coupling with Alfven modes with 
three extrema along field lines is allowed, we find no signatures of this resonant coupling in the 
power spectrum analysis nor the wave energy density evolution. 



5.2. Propagation of coupled fast and Alfven disturbances in a potential coronal 

arcade 

Oscillations in coronal magnetic structures are believed to be generated by nearby disturbances, 
such as flares or filament eruptions. It is clear that such disturbances are far from being a normal 
mode of a particular structure as our potential arcade. Therefore, we have next considered the 
impulsive excitation of perturbations by means of a localized disturbance, which is expected to be 
a better representation of the real phenomena that often trigger waves and oscillations in the solar 
corona. In particular a Gaussian velocity perturbation is considered and the response of the system 



is expected to be different from the one described in § 5.1 since now the initial perturbation is 
likely to be decomposed in a linear combination of normal modes with different frequencies that 
will constitute the resulting propagating wave. 

We have produced an impulsive excitation of the v \ n velocity component of the form given by 



Equation ( 15 ) and have considered k y ^ 0. Time evolution of the velocity components and the total 



energy density are displayed in Figure [7j which shows that the generated wave has both normal 



and perpendicular velocity components. Note that v\ y = in the absence of k y (see Terradas 



et al. 2008). It is clear in Figure [7] that the perturbed normal velocity component evolution is 



similar to the one presented by Terradas et al. (2008), for the decreasing Alfven speed model with 
constant density and k y = 0. For the normal velocity component, the shape of the wavefront is 
not circular, due to the fact that perturbations propagate faster toward the photosphere. For large 
times the front tends to be planar as the initial curvature of the wave packet is lost. As for the 
perpendicular velocity perturbation that is excited because of the three-dimensional character of 
the wave, its spatial distribution is highly anisotropic, with the signal concentrated around many 
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magnetic surfaces. A wavefront with fast-like properties, similar to the one present in vi n , can also 
be seen to propagate upwards producing v\ y perturbations until it leaves the system. At the end, 
a collection of Alfvenic oscillations are generated in the arcade. By comparing with the results 
presented in the previous section we can think about them as being generated by the resonant 
coupling between the fast- like wavefront and several Alfven continuum solutions, instead of the 
single resonance case shown in § |5.1| Once excited, magnetic surfaces remain oscillating with their 
natural period and for large times they are phase-mixed because of the transverse non-uniformity. 

We have next analyzed in a quantitative way the effect of k y on the properties of the gener- 
ated fast-like wavefront and the induced Alfvenic oscillations. Regarding the fast-like wavefront, 
Figure [8] (top-panels) shows different snapshots of the cut along x = of the v\ n component for 
different values of k y . These figures indicate that the larger the value of k y the faster the wavefront 
propagates. The propagation velocity can be measured by plotting the position of the wavefront 
maximum as a function of time (see Figure [8] bottom-left). The time evolution of the wavefront is 
followed for t ^ to and the initial position of the maximum is denoted by zq. For this relatively 
simple case, the numerical results can be compared with the analytical formula obtained by the 
integration of the local Alfven speed profile (see Equation [5] in Oliver et al. 1998). The resulting 
expression is 
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where the + and — signs correspond to upward and downward propagation, respectively. Figure [8] 
(bottom-panels) shows a perfect correspondence between the numerically measured speed and the 
analytical expression when different models of the solar atmosphere are considered. The increase 
of the travel speed of fast-like wavefronts when k y ^ is an important property to be taken into 
account in the three-dimensional problem. 

For the Alfvenic oscillations the power spectrum is analyzed in a cut along x = 0, which allows 
us to study the power on different magnetic surfaces. Figure [9] shows power at a large number of 
magnetic surfaces, not just around a selected group of field lines around a given magnetic surface, 
so a wide range of magnetic surfaces are excited because of the coupling. Also, not just the 
fundamental mode is excited, but also several higher harmonics. All of them have even parity with 
respect to x = 0, such as corresponds to the parity of the v\ n perturbation and the parity rule for 



k y ^ (Arregui et al. 2004a). When comparing the power spectrum obtained from the numerical 



solution with the analytical Alfven continuum frequencies given by Oliver et al. (1993) for the case 
k y = 0, we see that the signal coincides with the analytical curves for k y = 0. Such as expected, 
perpendicular propagation has no effect on the frequencies of Alfven waves generated on different 
magnetic surfaces in the arcade. This is a known result since oj ~ (k • B). 

As with the normal mode case, a quantitative analysis of the time-evolution of the wave energy 
of the system helps to better understand the process of energy conversion between fast and Alfven 
waves. Figure [7] (right-hand side panels), shows the evolution of the total energy density as a 
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function of time. At early stages this quantity shows a clear signature of a fast-like wavefront 
propagating through the domain. For long times the energy deposition is spatially distributed on 
the whole system, not only on a single magnetic surface, as in the previous section. Although a 
large part of the energy leaves the system in the form of fast-like wavefronts, part of the energy 
remains trapped in the Alfvenic oscillations that are resonantly excited in the arcade. The amount 
of energy trapped in the system can be calculated by the integration of the total energy density (see 
Equation [18]) in the whole domain as a function of time. The result is shown in Figure [To] (solid 
line). For short times the total energy remains almost constant, but when the fast front reaches 
the boundaries of the system a strong decrease of this quantity is seen. Resonant wave conversion 
very quickly produces velocity perturbations in the y-direction and the energy associated to these 
Alfvenic components grows up to its maximum value before the fast wavefront leaves the system. 
At later stages, a fraction of around 5 — 10% of the initial total energy is retained in the system 
and the total energy remains almost constant in the subsequent time evolution. We must note that 
this energy is trapped even in the absence of any density enhancement or wave cavity. 

So far, we have used fixed values for the perpendicular propagation wavenumber, k y , and 
the ratio of magnetic to density scale heights, 5. We have next analyzed the influence of these 
parameters on the obtained results concerning the energy transfer between fast and Alfven waves. 



Figure 11 shows the total energy time evolution for different values of k y . Several conclusions can 
be extracted. First, the amount of energy that is trapped by the system in the form of Alfvenic 
oscillations increases with k y and is above 40% for the largest value of this parameter that we have 
considered. This can be understood in terms of stronger resonant coupling occurring for larger 
values of k y . The relation between the total energy and the energy associated to the y-direction 
also changes with k y , in such a way that, while for relatively small k y almost all the energy of the 
system is stored in oscillations in the y-direction, for larger values of k y there is a difference between 
the total energy and the Alfvenic energy for large times. To understand this we need to mention 
that for k y ^ Alfven waves have both perpendicular and normal velocity components and so 
Alfven wave energy is not only contained in the y-direction. Although this effect is less visible in 



the simulations it can be measured, such as shown in Figure 11 Note also that for large times the 
two energy densities decay. This is due to numerical damping, since when very small scales are 
created the spatial resolution used is not fine enough to handle the localized Alfvenic oscillations 



that are phase- mixed for large times (see § 4.2.2) 



As for the 5 parameter, it controls our model atmosphere, since it allows us to select different 
ratios of the magnetic scale height to the density scale height. By repeating the previous numerical 
experiments for two additional values of this parameter the following results are obtained (see 
Figure [8] bottom-panels). Depending on the value of 5, the Alfven speed profile in the vertical 
direction has a steeper or flatter profile. This means that the time that a fast-like perturbation 
needs to reach the boundaries of the system and leave it varies with 5. Therefore, the time at 
which the sudden decrease of the total energy of the system occurs differs for different values of 5, 



see Figure 10 However, the fractional amount of wave energy that is transferred to Alfven waves 
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and is trapped in the system does not depend on the model atmosphere we use. Nevertheless, the 
rate at which the energy transfer occurs does depend in the model atmosphere, such as can be 



appreciated from the different slopes of the energy in Figure 10 



6. Conclusions 

In this paper we have studied the temporal evolution of coupled fast and Alfven waves in a 
potential coronal arcade when three-dimensional propagation is allowed. Because of the inclusion 
of three-dimensional dependence on the perturbed quantities, fast and Alfven waves are coupled 
and the resulting solutions display a mixed fast/ Alfven character. The non-uniform nature of the 
considered medium produces the coupling to be of resonant nature, in such a way that transfer of 
energy and wave damping occur in the system. 

First, the nature of resonant coupling between a fast normal mode of the system and Alfven 
continuum modes has been analyzed. It is seen that the fast mode with a global nature resonantly 
couples to localized Alfven waves around a given magnetic surface in the arcade, thus transferring 
its energy to the later. The position of the resonant surface perfectly agrees with the resonant 
frequency condition predicted by several authors in previous studies of this kind, and with the 



parity rules given by Arregui et al. (2004a). 



Next, the temporal evolution of a localized impulsive disturbance has been analyzed. The 
inclusion of perpendicular propagation produces an increase in the wave propagation speed for the 
fast-like wavefront when compared to the purely poloidal propagation case. As in the previous case, 
perpendicular propagation induces the excitation of Alfvenic oscillations around magnetic surfaces, 
due to the resonant coupling between fast and Alfven waves. Now these oscillations cover almost 
the whole domain in the arcade, so that the energy of the initial perturbation is spread into localized 
Alfvenic waves. The frequency of the induced Alfvenic oscillations is seen to be independent from 
the perpendicular wavenumber. As time progresses and the initial wavefront leaves the system 
part of the energy is stored in these Alfven waves which remain confined around magnetic surfaces. 
Phase mixing then gives rise to smaller and smaller spatial scales, until the numerical code is unable 
to properly follow the subsequent time evolution. The energy trapping around magnetic surfaces 
occurs even in the absence of a density enhancement or a wave cavity structure, and is only due to 
the non- uniformity of the density profile and the magnetic structuring, which lead to a non-uniform 
Alfven speed distribution. 

Finally, the efficiency of this wave energy transfer between large scale disturbances and small 
scale oscillations has been studied as a function of the perpendicular wavenumber and for different 
values of the ratio of the magnetic scale height to the density scale height. It is seen that the 
first factor strongly affects the amount of energy trapped by Alfven waves. The amount of energy 
trapped by the arcade increases for increasing value of the perpendicular wavenumber. The partic- 
ular ratio of magnetic to density scale heights determines how fast the available fast wave energy 
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leaves the system and, therefore, the rate at which energy can be transferred to Alfven waves, but 
not the final amount of energy stored by the arcade in the form of Alfvenic oscillations. These 2.5D 
simulations should be extended in the future to more realistic 3D simulations in order to ascertain 
the applicability of our conclusions to the real wave dynamics observed in coronal structures. 
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acknowledges MCyT for a fellowship. 

REFERENCES 

Arregui, I., Oliver, R., & Ballester, J. L. 2004a, A&A, 425, 729 
-. 2004b, ApJ, 602, 1006 

Aschwanden, M. J., de Pontieu, B., Schrijver, C. J., & Title, A. M. 2002, Sol. Phys., 206, 99 

Aschwanden, M. J., Fletcher, L., Schrijver, C. J., & Alexander, D. 1999, ApJ, 520, 880 

Banerjee, D., Erdelyi, R., Oliver, R., & O'Shea, E. 2007, Sol. Phys., 246, 3 

Bona, C, Bona-Casas, C, & Terradas, J. 2009, Journal of Computational Physics, 228, 2266 

Brady, C. S., & Arber, T. D. 2005, A&A, 438, 733 

Brady, C. S., Verwichte, E., & Arber, T. D. 2006, A&A, 449, 389 

Cadez , V. M., & Ballester, J. L. 1995a, A&A, 296, 537 

-. 1995b, A&A, 296, 550 

Del Zanna, L., Schaekens, E., & Velli, M. 2005, A&A, 431, 1095 
Diaz, A. J., Zaqarashvili, T., & Roberts, B. 2006, A&A, 455, 709 
Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 
Goossens, M., Poedts, S., & Hermans, D. 1985, Sol. Phys., 102, 51 
Halberstadt, G., & Goedbloed, J. P. 1993, A&A, 280, 647 
Murawski, K., Selwa, M., & Nocera, L. 2005, A&A, 437, 687 
Nakariakov, V. M., & Ofman, L. 2001, A&A, 372, L53 

Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., & Davila, J. M. 1999, Science, 285, 862 



-16- 



Oliver, R., Ballester, J. L., Hood, A. W., & Priest, E. R. 1993, A&A, 273, 647 

Oliver, R., Hood, A. W., & Priest, E. R. 1996, ApJ, 461, 424 

Oliver, R., Murawski, K., & Ballester, J. L. 1998, A&A, 330, 726 

Poedts, S., & Goossens, M. 1988, A&A, 198, 331 

Poedts, S., Hermans, D., & Goossens, M. 1985, A&A, 151, 16 

Roberts, B., Edwin, P. M., & Benz, A. O. 1984, ApJ, 279, 857 

Ruderman, M. S., Goossens, M., Ballester, J. L., & Oliver, R. 1997, A&A, 328, 361 

Russell, A. J. B., & Wright, A. N. 2009, in press 

Schrijver, C. J., Aschwanden, M. J., & Title, A. M. 2002, Sol. Phys., 206, 69 

Selwa, M., Murawski, K., Solanki, S. K., & Wang, T. J. 2007, A&A, 462, 1127 

Selwa, M., Solanki, S. K., Murawski, K., Wang, T. J., & Shumlak, U. 2006, A&A, 454, 653 

Terradas, J. 2009, Space Science Reviews, in press. doi:10.1007/sll214-009-9560-3 

Terradas, J., Oliver, R., & Ballester, J. L. 1999, ApJ, 517, 488 

-. 2006, ApJ, 650, L91 

Terradas, J., Oliver, R., Ballester, J. L., & Keppens, R. 2008, ApJ, 675, 875 

Tsubaki, T. 1988, in Solar and Stellar Coronal Structure and Dynamics, ed. R. C. Altrock, 140-149 

Van Doorsselaere, T., Debosscher, A., Andries, J., & Poedts, S. 2004, A&A, 424, 1065 

Van Doorsselaere, T., Verwichte, E., & Terradas, J. 2009, Space Science Reviews, in press. 
doi:10.1007/sll214-009-9530-9 

Verwichte, E., Foullon, C, & Nakariakov, V. M. 2006a, A&A, 446, 1139 

-. 2006b, A&A, 449, 769 

-. 2006c, A&A, 452, 615 

Verwichte, E., Nakariakov, V. M., Ofman, L., & Deluca, E. E. 2004, Sol. Phys., 223, 77 
Wright, A. N. 1992, J. Geophys. Res., 97, 6429 



This preprint was prepared with the AAS MgX macros v5.2. 



-17- 




Fig. 1. — Sketch of the magnetostatic configuration of a potential coronal arcade, where the solid 
curves represent magnetic field lines, given by A(x,z) = constant. These curves in the xz-plane 
become arcade surfaces in three dimensions. In this model z measures the upward distance from 
the base of the corona (placed at z = 0). The three orthogonal unit vectors, e n , ej_ and ey, defining 
the normal, perpendicular and parallel directions respectively, are also shown at a particular point. 
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Fig. 2. — (a) Shaded contours are the power spectrum of v\ y at x = as a function of height, z/L, 
and normalized frequency, uL/vao- In this simulation the spatial grid is set to 600 x 600 while 
the numerical dissipation is fixed to a n = 0.001. Solid lines are the theoretical frequency of the 
normal Alfven modes for 5 = 0, given by Oliver et al. (1993). From bottom to top they represent 
the fundamental mode and its harmonics, (b) Temporal evolution of the v\ y component when the 
initial perturbation is an Alfven normal mode of the system located at x = and z m /L = 0.33 (see 
Equation [l7]). Here 5 = 0, the numerical grid has 400 x 400 points, and a n = 0.001. Different 
solutions correspond to a/L = 0.01 (solid), a/L = 0.1 (dash-dotted) and a/L = 0.2 (long-dashed). 
Time is given in units of ta = L/vao- 
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Fig. 3. — (a) Temporal evolution of the v\ y velocity component at x = and z m /L = 0.33 for 



an initial perturbation given by Equation (17) with a/L = 0.2 and a n = 0.001. (b) Temporal 
evolution of the v\ y velocity component at x = and z m /L = 0.66 for an initial perturbation with 
a/L = 0.2 and a n = 0.001. In both panels solid, dash-dotted, and long-dashed lines represent a 
resolution of 200 x 200, 400 x 400, and 600 x 600 points, respectively. 
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Fig. 4. — Several snapshots of the two velocity components, v\ n and v\ y , and of the total energy 
density in a potential arcade with 5 = 2 and for k y L = 1. The initial perturbation in v\ n is 
the fundamental normal mode of the k y = problem. For this simulation a 600 x 600 grid is 
used. Magnetic field lines are represented with white (left and middle panels) and black lines (right 
panels). A movie displaying the full time evolution is available in the electronic version of the 
journal. 
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Fig. 5. — Temporal evolution of the normal, v\ n (solid line), and perpendicular, v\ y (dashed line), 
velocity components at x = 0, z/L = 0.35. Data taken from the simulation shown in Figure Hi 
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Fig. 6. — The shaded contours represent the power spectrum of (a) 5E n and (b) SE y as a function 
height, z/L, at the symmetry plane x = 0, for the simulation shown in Figure |4j Note that because 
of the quadratic nature of the wave energy density the curved lines have double the frequency of 
the Alfven continua given by Oliver et al. ( 1993 ) and the horizontal lines have double the frequency 
of the fast normal mode. 
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Fig. 7. — Time evolution of the v\ n (left panels), v\ y (middle panels) velocity components, and of 
the total energy density (right panels) in a potential arcade with 5 = 1. The initial perturbation 



is imposed on the v\ n component with x s = 0, z s = 1, v s = lO -4 ?;^) o, = 0.2L (see Equation 15) 
and k y L = 1. For this simulation a 600 x 600 grid is used. Magnetic field lines are represented 
with white (left and middle panels) and black lines (right panels). A movie displaying the full time 
evolution is available in the electronic version of the journal. 
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Fig. 8. — Top-panels Several snapshots of the spatial distribution of the normal velocity compo- 
nent along x = 0. The values of the longitudinal wavenumber are k y = (solid line), k y L = 3 
(dotted), k y L = 7 (dashed), and k y L = 10 (dash-dotted). Bottom-panels Position of the wave- 
front as a function of time when different values of the delta parameter, 5 = 1 (left), 6 = 2 (middle), 
and 5 = 3 (right), and longitudinal wavenumber, k y = (squares) and k y L = 10 (triangles), are 
selected. The solid line shows the analytical solution for the wavefront position when longitudinal 
propagation is not allowed (k y = 0); see Equation (19). 
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Fig. 9. — Shaded contours represent the normalized power spectrum of the v\ y velocity component 
corresponding to the simulation show in Figure [7] as a function of the maximum height of field 
lines, z/L, and normalized frequency, uL/vao- Solid lines are the theoretical frequency of the 
Alfven normal mode obtained by Oliver et al. (1993). The frequency analysis is made at the 
symmetry plane, x = 0. 
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Fig. 10. — The solid (5 = 1), dotted (5 = 2), and dashed (5 = 3) lines are the normalized total 
energy (thick lines) and the normalized total energy associated to the y-direction (thin lines) as a 
function of time for k y L = 1. 
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Fig. 11. — Normalized total energy of the system as a function of time (solid line) and normalized 
energy associated to the y-direction (dotted line) for 5 = 1 and different values of the longitudinal 
wavenumber. (a) k y L = 3, (b) k y L = 5, (c) k y L = 7, and (d) k y L = 9. 



